A search for high-mass resonances decaying to dimuons at CDF 
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We present a search for high-mass neutral resonances using dimuon data corresponding to an 
integrated luminosity of 2.3 fb^^ collected in pp collisions at = 1-96 ToV by the CDF II detector 
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the dimuon invariant-mass spectrum. We set 95% confidence level upper limits on a-BR{pp ^ X — > 
fj,jj,), where X is a boson with spin 0, 1, or 2. Using these cross section limits, we determine lower 
mass limits on sneutrinos in Ji-parity- violating supersymmetric models, Z' bosons, and Kaluza-Klein 
gravitons in the Randall-Sundrum model. 
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Neutral resonances decaying to muons have historically 
been a source of major discoveries. They also occur in a 
variety of theoretical models which attempt to unify the 
standard model (SM) forces or explain the large gap be- 
tween the SM and gravitational energy scales. The gauge 
group SU(3)c X SUi2)L x C/(l)y of the SM can be em- 
bedded in larger gauge groups such as SU{5), 50(10), 
and Eq, to achieve unification in a grand unified theory 
(GUT) 0,0,13,11. In many schemes of GUT symmetry- 
breaking, C/(l) gauge groups survive to relatively low en- 
ergies (4I , leading to the prediction of neutral gauge vec- 
tor {Z') bosons. Such Z' bosons typically couple with 
electroweak strength to SM fermions, and can be ob- 
served at hadron colliders as narrow, spin-1, dimuon res- 
onances from qq ^ Z' /ifl. Many other models, such 
as the SU{2)l x SU{2)r x U{\)b-l gauge group of the 
left-right model [i], and the "little Higgs" models [1,0], 
also predict heavy neutral gauge bosons. 

Additional spatial dimensions are a possible expla- 
nation for the gap between the electroweak symmetry- 
breaking scale and the gravitational energy scale 
Mpianck B il- The Randall-Sundrum (RS) model Q 
predicts excited Kaluza-Klein modes of the graviton, 
which appear as spin-2 resonances G* in the process 
qq G* fip,. These modes have a narrow intrinsic 
width when fc/Mpianck < 0.1, where is the spacetime 
curvature in the extra dimension. In superstring theories 
with 0(1) couplings, fc/Afpia„ck ~ 0.01 [l3|. 

Spin-0 resonances such as the sneutrino i> in the pro- 
cess gg — > z/, — > /i/I are predicted by supersymmetric 
theories with i?-parity violation [Hj . Scalar Higgs bosons 
can be produced as resonances and decay to dimuons. 

The most sensitive direct searches for high-mass bo- 
son resonances, which have previously been performed at 
the Tevatron, have set 95% confidence level (C.L.) lower 
limits on the masses Mz', Mq*, and Mp of Z' bosons, 
RS gravitons, and sneutrinos, respectively. The previous 
dimuon pubhcation from CDF II, based on ^ 200 pb"^ 
of integrated luminosity [13, set mass limits that vary 
from 170 GeV to 885 GeV[li| depending on the boson 
spin and couplings to the SM fermions. Other dilepton 



and diphoton decay channels have also been explored at 
the Tevatron IJ, [15| . Using an order of magnitude more 
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data, we present in this Letter the most sensitive direct 
search to date for Z', G*, and i> bosons at high mass. 
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FIG. 1: The distribution of (TeV"^) for the observed 

data (points), the individual backgrounds (dotted or dashed 
histograms) and the summed background (solid histogram). 
The Z boson peak is prominently seen. The inverse mass 
distribution has the useful feature that the detector resolution 
is constant (« 0.17 TeV~^) over the range shown in the plot. 



This analysis uses 2.3 fb ^ of data fro m pp collisions at 
^ = 1.96 TeV in the CDF II detector [ll[l3]- CDF II 
is a magnetic spectrometer surrounded by calorimeters 
and muon detectors. We use the central drift cham- 
ber (COT) the central calorimeter [l^, and the 
muon detectors for identification and measurement 
of muons with jryj < 1 [l3| . The online selection requires 
a COT track with px > 18 GeV [l^, and matching muon 
detector hits. We select a pair of oppositely-charged 
muons, each with a COT track with px > 30 GeV passing 
quality requirements, and a minimum-ionization signal in 
the calorimeter. Cosmic rays are rejected using COT hit 
timing ^)\. The dimuon signal sample consists of 68150 
events in the control region 70 < m^p < 100 GeV, where 
the pp ^ Z — > fip, process dominates, and 3804 events in 
the search region m^p > 100 GeV. 

The alignment of the COT is performed using a pure 
sample of high-momentum cosmic-ray muons, in order to 
obtain the best possible dimuon mass resolution. Each 
muon's complete trajectory is fitted to a single helix [2l[. 
The fits are used to determine the relative locations 
of the sense wires, including gravitational and electro- 
static displacements, with a statistical accuracy of a few 
microns |17| . We constrain remaining misalignments, 
which cause a bias in the track curvature, by compar- 
ing (E/p) [l3| for electrons and positrons. The tracker 
momentum scale and resolution is measured by template- 
fitting the Z —f fj,p, mass peak, and calibrating to the 
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world average values (24] of the Z boson mass and width. 
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FIG. 2: The difFerence between the distributions of ra,,r, 
(TeV~^) for the observed data and the summed background, 
divided by the expected statistical uncertainty in each bin. 
All vertical error bars have unit size. 

For a resonance with electroweak coupling and mass 
above 200 GeV, the observed width of the m^^ dis- 
tribution is dominated by the track curvature resolu- 
tion, resulting in an approximately constant resolution 
0.17 TeV~^. Our search strategy is to con- 



struct templates of the observable Tn~^ distribution for a 
range of boson Breit-Wigner pole masses, add the back- 
ground distributions to the templates, and compare the 
templates to the distribution from the data in the 
search region m^^ > 100 GeV. The simulated templates 
(including backgrounds) are normalized to the data in 
the 70 GeV < m^p < 100 GeV region, thus cancelhng 
several sources of systematic uncertainty. 

We determine the most likely number of signal events 
(A^s)j and the corresponding confidence intervals [23t . 
from the binned Poisson likelihood [l3| for the observed 
data to be produced by a sum of signal and background 
templates. The use of the constant-resolution variable 
TO~p simplifies the optimization of the template binning 
and the scan over the boson pole masses. 

Signal and SM Drell-Yan background distributions are 
evaluated using a specialized Monte Carlo (MC) simu- 
lation [l3| of boson production and decay, and of the 
detector response to the leptons and hadrons. The kine- 
matics of boson production and decay are obtained from 
the PYTHIA [il] event generator using the CTEQ6M [25| 
set of parton distribution functions. QED radiation is 
simulated [l3] based on the WGRAD program The 
MC performs a detailed hit-level simulation of the lep- 
ton tracks. COT hits are generated according to their 
resolution 150 fxro) and measured efficiencies, and a 
helix fit is performed (as it is in data) to simulate the re- 
constructed track. We apply a mass-dependent next-to- 
next-to-leading order (NNLO) multiplicative correction 
(if-factor) [27| to the SM Drell-Yan background. 

The SM production processes for W^W~ [2^ and 
tt [2^ have small contributions, and are evaluated using 
their NLO cross sections, PYTHIA, and a detector simula- 
tion based on GEANT [30| . Misidentification backgrounds 
result from cosmic rays, QCD jets, and -k/K decays-in- 
flight (DIF). We evaluate the cosmic-ray background us- 
ing a large sample of cosmic rays identified with the COT- 



timing-based algorithm [2l|, and using the direction-of- 
fiight information provided by this algorithm. The ra~^ 
shape of misidentified jets is evaluated from a large sam- 
ple of inclusive jet events. Decays-in-fiight within the 
COT active volume generate a kink along the helical tra- 
jectory, resulting in a mismeasurement of the track cur- 
vature. For large reconstructed momenta, the measured 
DIF curvature distribution is approximately uniform and 
m-;(Tev')io leads to a flat m~p spectrum. Most DIF tracks are re- 
jected using their abnormal COT-hit pattern and large fit 
X^. The jet and DIF backgrounds are normalized using 
the mass distribution of same-charge dimuon events. 

Figure [T] shows the distributions of the observed 
data and the expected backgrounds, which are in good 
agreement (as shown in Fig. [2|) . A resonance whose ob- 
served width is dominated by detector resolution would 
appear as a peak spanning approximately three bins. 
The likelihood-based fitter finds no significant excess. 
We use background-only ensembles of simulated events, 
each with the statistics of the data sample, to evalu- 
ate the probability of statistical fluctuations anywhere 
in the search region generating a discrepancy at least as 
significant as the largest discrepancy found in the data. 
We find this probabihty ("p- value") to be 6.6% and we 
conclude that the observed data are statistically consis- 
tent with the SM expectation. The dielectron niee spec- 
trum from 2.5 fb~^ of CDF II data [3l| shows that the 
largest discrepancy with the expected background occurs 
at TTiee ~ 240 GeV. Figure [2] shows that the dimuon data 
are consistent with the expectation near this mass to bet- 
ter than Ict in statistical precision. The sensitivity of the 
dielectron analysis for a spin-1 resonance at this mass is 
w 20% better than the dimuon analysis reported here. 

The likelihood fitter determines the 95% C.L. upper 
limit on the number of signal events, for each value of 
the resonance pole mass. We convert these limits to 
limits on a ■ BR(y,D — > /i/2), <J ■ BR{Z' — > /i/l), and 
tJ • BR{G* jijl) using the total acceptance as a func- 
tion of pole mass, the NNLO cross section for Z ^fi 
of 251.3 pb [iBl, and dividing by the observed number of 
Z ^ ^p, events. The acceptance is verified with the de- 
tailed GEANT-based simulation, and comparisons to data 
distributions. The muon identification efficiency is veri- 
fied using a pure data sample of Z bosons triggered by 
one identified muon. The total acceptance, including 
kinematic and fiducial acceptance and dimuon identifi- 
cation, increases from « 13% (« 20%) for a pole mass of 
90 GeV to « 40% (« 45%) for a Z' (graviton) pole mass 
of 1 TeV, and decreases for higher pole masses due to the 
kinematic limit of the parton collisions. The 95% C.L. 
upper limits on cr • BR{i>,i? fJ-p.), o ■ BR{Z' /i/l), 
and cr • BRiG* fip,) are shown in Fig. [3] The domi- 
nant mass-dependent systematic uncertainties arise from 
parton distribution functions (16%), the NNLO if-factor 
(9%) [131, QED radiative corrections (3%) and ac- 
ceptance (3%), all quoted at 1 TeV. These uncertainties 
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FIG. 3: The 95% C.L. upper limits on a ■ BR{i>, i -> plJj.) 
vs Mc (top), a ■ BR{Z' /ijl) versus Mz' (middle), and 
a ■ BR{G* fip,) versus Mc (bottom). Also shown are the 
theoretical cross sections for various model parameter val- 
ues [1, [ll|, [13 • The expected limits and ranges of limits, 
as derived from simulated experiments (SE), are shown for 
comparison. The step size between adjacent templates in the 
signal scan is 0.2 TeV~^ in pole mass. 

are incorporated as functions of m^p and increase mono- 
tonically beyond 100 GeV. Uncertainties on the momen- 
tum scale and resolution, and on the non-Drell-Yan back- 
ground predictions, have a negligible effect. 

Our signal templates have been generated with a res- 
onance pole width T = 2.8% x M, based on the SM Z 
boson width. Thus our signal scan probes an observed 
width of « [17%(M/TeV) ® 2.8%] M. In a model where 
the observed width increases by a factor x, the cross sec- 
tion limits would increase by about a factor of ^/x. 

We use PYTHIA to compute the cross sections for pro- 
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RS graviton 


graviton 
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model 


mass limit 


k/ i\/pianck 
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Z'j 


789 


0.01 


293 


0.0001 


397 


7' 


821 


0.015 


409 


0.0002 


441 




861 


0.025 


493 


0.0005 


541 


Zip 


878 


0.035 


651 


0.001 


662 




892 


0.05 


746 


0.002 


731 


Z^ 


904 


0.07 


824 


0.005 


810 




1030 


0.1 


921 


0.01 


866 



TABLE I: 95% C.L. lower limits on Z' , graviton, and sneu- 
trino masses (in GeV) for various model parameters [9l.[Tll.[33. 
For the _R-parity-violating sneutrino model, A is the ddii cou- 
phng and BR denotes the u,v jift branching ratio. 



duction of Z' bosons predicted by Eq models [33j or hav- 
ing the same couplings to SM fermions as the Z boson, 
and of G* bosons for various fc/Afpianck values. We ap- 
ply the NNLO if-factor to these LO cross sections. The 
NLO z> production cross sections are obtained from 
We derive the boson mass limits shown in Table HI 

In conclusion, we have presented a direct search for 
high-mass neutral resonances with spin-0, 1, and 2, us- 
ing an integrated luminosity of 2.3 fb~^ collected by the 
CDF II detector. Our dimuon invariant mass spectrum is 
consistent with the SM expectation. We set the world's 
tightest constraints on Z' bosons in various models, on 
Kaluza-Klein graviton modes in the RS model, and on 
sneutrinos in i?-parity violating supersymmetric models. 
At 95% C.L., we exclude 100 < Mz' < 982 GeV for 
a Z'^ boson of the Eg model, 100 < Mg- < 921 GeV 
for fc/Mpia„ck = 0.1, and 100 < Mc < 810 GeV for 
• BR{v, V ^p.) = 0.01, where A is the ddi> coupHng 
and BR denotes the v,!? ^ fj,fl branching ratio. 
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